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A common approach for the numerical simulation of wave propagation on a spatially unbounded 
domain is to truncate the domain via an artificial boundary, thus forming a finite computational 
domain with an outer boundary. Absorbing boundary conditions must then be specified at the 
boundary such that the resulting initial-boundary value problem is well posed and such that the 
amount of spurious reflection is minimized. In this article, we review recent results on the con- 
struction of absorbing boundary conditions in General Relativity and their application to numerical 
relativity. 

o 

O ■ I. INTRODUCTION 

To numerically simulate the evolution of hyperbolic partial differential equations on a spatially unbounded domain, 
i one usually replaces the unbounded domain with a finite, compact domain Vl with artificial outer boundary <90. 
Boundary conditions on 9f2 must then be specified in order to obtain a unique Cauchy evolution. These conditions 

■ should be formulated so that they form a well posed initial boundary value problem (IBVP) and, ideally, are completely 
£f) ■ transparent to the physical problem on the unbounded domain, fn practice, complete transparency cannot be achieved 

easily and therefore, the boundaries introduce spurious reflections into the computational domain CI. The idea then is 
to develop what is called absorbing boundary conditions which form a well posed IBVP and insure that the amount 
qt< of spurious reflection from dQ is as small as possible. 

There has been a substantial amount of work on the construction of absorbing (also called non-reflecting in the 
j^j^' literature) boundary conditions for wave problems in acoustics, electromagnetism, meteorology, and solid geophysics 
(see [l| for a review). One approach is based on a hierarchy of local boundary conditions 0, H, 0] with increasing 
order of accuracy. Although higher order local boundary conditions usually involve solving a high order differential 
equation at the boundary, the problem can be dealt with by introducing auxiliary variables at the boundary surface 
Q . A different approach is based on fast converging series expansions of exact nonlocal boundary conditions (see 
and references therein). 

Constructing absorbing outer boundary conditions in General Relativity is difficult. First of all, Einstein's field 
equations determine the evolution of the metric tensor, so one does not know the geometrical structure of the spacetime 
before actually solving the IBVP. Hence, it is not clear a priori how the geometry of the outer boundary evolves. 
One could turn the argument around to guide the choice of boundary conditions so that they fix the embedding of 
C""^ , the boundary surface into the resulting spacetime. Second, in the Cauchy formulation of Einstein's field equations, 

■ one usually encounters constraint fields which propagate across the boundary into the computational domain. This 
is in contrast to the standard Cauchy formulation of Maxwell's equations, where the evolution equations imply that 
the constraint variables (namely, the divergence of the electric and magnetic fields) are constant in time. As a 
consequence, in the Einstein case, constraint-preserving boundary conditions need to be specified at dfl such that 
constraint violations are not introduced into the computational domain. Finally, due to the nonlinear nature of the 
theory and its diffeomorphism invariance, it is difficult to define precisely what is meant by in- and outgoing radiation 
in General Relativity. Therefore, it is not even clear how to quantify the amount of spurious gravitational reflection 
from the boundary. These issues all contribute to the challenge of developing accurate absorbing boundary conditions 
for Einstein's field equations. 

In this article, we briefly review some recent results on the construction of absorbing outer boundaries in General 
Relativity. We start with a discussion of absorbing outer boundaries for the wave equation on Minkowski and 
Schwarzschild spacetimes in Sect. |TI]where we review results from the literature and derive asymptotic outgoing wave 
solutions on weakly curved spacetimes. Recent attempts to generalize these results to Einstein's field equations are 
mentioned in Sect. IIII1 In particular, we discuss constraint-preserving boundary conditions, boundary conditions 
designed to minimize spurious reflections of gravitational waves, recent results on well posed initial-boundary value 
formulations and recent applications to numerical relativity. Concluding remarks are drawn in Sect. IIV1 

There exist other approaches for dealing with gravitational wave propagation on an infinite domain which are not 
discussed in this article. These involve matching techniques (see Ref. B for a review) or avoiding introducing an 
artificial outer boundary altogether by compactifying spatial infinity 0, [2], or making use of hyperboloidal slices and 
compactifying null infinity (see, for instance, [Til Il2l Il3j). 
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II. ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION 



To illustrate some of the ideas involved in constructing absorbing boundary conditions, we start with three examples 
in order of increasing difficulty. The simplest example is the flat wave equation on an interval (—1, 1). In this case, 
perfectly absorbing boundary conditions can be constructed since the most general solution is simply the superposition 
of an arbitrary function of retarded time with an arbitrary function of advanced time. The second example discusses 
the flat wave equation on a three-dimensional ball Br of radius R > 0. This example is already much more difficult 
than the previous one since waves can now travel in infinitely many directions. One strategy here is to obtain a 
hierarchy of absorbing boundary conditions with increasing order of accuracy. The last example generalizes the 
second example to a weakly curved background. When the curvature of the background is taken into account, a 
wave solution propagating on Minkowski spacetime acquires two types of correction terms. The first is a curvature 
correction term which obeys Huygens' principle and the second is a fast decaying term which violates Huygens' 
principle and describes the backscatter off curvature. It is important to understand these effects before constructing 
absorbing boundaries. An accurate boundary condition must not eliminate the backscatter which is a real, physical 
effect! This example shows that an absorbing boundary condition (as defined in the introduction) does not mean that 
all reflections off the boundary are to be eliminated, but only spurious reflections. 

A. The flat wave equation on an interval 

Consider the one-dimensional flat wave equation 

(d 2 t - d 2 x ) u(t, x) = 0, t>0, x€ [-1,1]. (1) 
The general solution is a superposition of a left- and a right-moving solution, 

u(t,x) = f^(x + t) + fXx-t). 

Therefore, the boundary conditions 

(d t -d x )u(t,-l) = 0, (d t +d x )u(t,+l) = 0, t>0, (2) 

are perfectly absorbing according to our terminology. Indeed, the operator b\ := dt + d x has as its kernel the right- 
moving solutions j ' /-{x — t), hence the boundary condition biu(t, 1) = Q, t > 0, does not "touch" these solutions. 
On the other hand, bif\^(t + x) = 2fl(t + x), so the boundary condition at x = 1 requires that f\(v) = f\(l) is 
constant for advanced time v = t + x > 1. A similar argument shows that the left boundary condition implies that 
f y(—u) — //•(— 1) is constant for retarded time u = t — x > 1. Furthermore, it is known that the conditions (J2|) 
together with Eq. (TT]) and suitable initial conditions for u and dtu at t — yield a well posed IBVP. In particular, 
the solution is identically zero after one crossing time (t > 2) if the initial data has compact support. 

B. The flat wave equation on a three-dimensional ball 

Next, consider the three-dimensional flat wave equation 

(d 2 t - d 2 x - d 2 y - 81) u(t, x)=0, t>0, (x, y, z) e B R , (3) 
on a ball Br of radius R > 0. The general solution can be decomposed into spherical harmonics 
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which yields the family of reduced equations 



u im {t,r)=0, t >0, r G (0,R). (4) 



For i = this equation reduces to the previous example, and the general solution is woo(A r ) — Uqo /"( r ~ t) + Uoo^(r+t) 
with Uoo / and C/oox ^ w0 arbitrary functions. Therefore, the boundary condition 

6(ru)| P=fl = 0, (5) 
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where b :— r 2 (d t + d r ) is perfectly absorbing for spherical waves. For I > 1, exact solutions can be generated from 
the solutions for £ = by applying suitable differential operators to u^{t,r). For this, we define the operators [3| 



ag = d r 



which satisfy the operator identities 



If Qy 



,2 , 



As a consequence, for each £ = 1,2, 3..., we have 



2 , ^+1) 



ala' 



<9 t 2 + aja^ 



d t + <4-i a <-i 



ft^ J ■ ■ ■ ft ^ 



= alal_ v ..al[d?-d 2 r ]. 



Therefore, we have the explicit in- and outgoing solutions 



u hn *^(t,r) = a\a\_ l ...a\V im {r + t) = ^(-1 

3=0 



(2t-j)\ 



Utr, 



(t,r) = ata]_ 1 ..4^ m (r-t) = ^(-l)^|^|(2r)^^«(r- i ) ) 



(6) 



(7) (7) 

where Vg m and L^ m are arbitrary smooth functions with j'th derivatives and U^i, respectively. In order to 



construct boundary conditions which are perfectly absorbing for ug n 
b = r 2 (dt + d r ) as above, then 



one first shows the following identity. Let 



b l+1 a\a\_ v ..a[U{r - t) = 



(7) 



for all I = 0, 1, 2, ... and all sufficiently smooth functions U. This identity follows easily from Eq. (J6j) and the fact that 



b i+1 (r k ) =k{k + l)--- (k + e)r k 
condition 



k+i+l 



if k S {0, -1, -2, -£}. Therefore, given L e {1, 2, 3, ...}, the boundary 
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fe L + 1 (ra)| 



= 0, 



(8) 



leaves the outgoing solutions with £ < L unaltered. Notice that this condition is local in the sense that its formulation 
does not require the decomposition of u into spherical harmonics. Furthermore, it was shown in 0] for domains 
which can be more general than Br that each condition Bl yields a well posed IB VP. By uniqueness this implies 
that initial data corresponding to a purely outgoing solution with I < L yields a purely outgoing solution (without 
reflections). In this sense, the condition Bl is perfectly absorbing for waves with £ < L. For waves with £ > L, one 
obtains spurious reflections; however, for monochromatic radiation with wave number k, the corresponding amplitude 
reflection coefficients can be calculated to decay as (fci?) -2 ( i+1 ) in the wave zone kR 3> 1 [lq . Furthermore, in most 
scenarios with smooth solutions, the amplitudes corresponding to the lower few €s will dominate over the ones with 
high £ so that reflections from high ^'s are unimportant. For a numerical implementation of the boundary condition 
B2 via spectral methods and a possible application to General Relativity see Ref. [l6j . 



C. The wave equation on a weakly curved background 

Next, we generalize the previous example to the wave equation on a weakly curved background. More specifically, 
we consider the wave equation on the far field region of an asymptotically flat spacetime. Such a spacetime has the 
form of Minkowski plus 1/r correction terms, where r is the areal radius. The monopolar correction is given by the 
M/r term of the Schwarzschild metric, where M represents the total (Arnowitt-Deser-Misner) mass of the spacetime. 
Therefore, to first approximation, we may assume that spacetime is described by the exterior of a Schwarzschild metric 
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of mass M. In outgoing Eddington-Finkelstein coordinates (t — r,r) the reduced wave equation on a Schwarzschild 
background is 



2M 
r 



(dt + d r ) 2 --(dt + d r ) + 



l + cr 



(9) 



Here, a is a parameter which in the present case is zero but is left arbitrary for future convenience. For M = 0, this 
equation reduces to the fiat reduced wave equation for which outgoing solutions have the form 



Hfi/ftr) = alal_ 1 ...a\U£ m (r - t) 

with a sufficiently smooth function Ut m . For the following, we also assume that Ug m is bounded and vanishes for 
sufficiently negative values of its argument. If M > 0, we seek an approximate outgoing solution in the far field region 
r « R > M of the form 



t„t „ti 



, s (o) / n 2M m , x ~ (2M\ 
u im y{t,r) = u\^{t,r) + —ulJ^it^r) + O I — I 



(10) 



Plugging this into Eq. ([9]) and expanding in M/R, we find that must satisfy 
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(d t + d r ) 2 --(3 f + d r ) + 



l + cr 



/ (0) 



3=0 



The solution to this equation can be written as the sum of two terms, 

The first term includes corrections from the curvature and obeys Huygens' principle. It has the form 

i 

4Z v H^) = R^(-iy^^(2ry-%ut 1) (r-t), 



(11) 



(12) 



3=0 



where the coefficients Cq, c%, q can be computed from the recursion relation 

2(1 + a) 



C( = 0, Q_l = 0, C(-2 



{£-!)£{£ +!){£ + 2) ' 



c.,_i = c,- 



2(^-i) (^-j) 2 + ^ 



j =£-2,^-3, ...1. 



(13) 



The second term is a fast decaying term which violates Huygens' principle and describes the backscatter off the 
curvature of the background. It has the form 



(backscatter) 



(2£V £ 2 + a ( 
(t, r) = R(l - St) ^ J e+1 U im (r -t) + R(l + S t a) J K e (t,r, x)U em (x)dx, (14) 



where St = 1 for £ = and St = for I > 1, and the integral kernel r, x) is given by 

1 1 ^(2£-j)l 



Kg(t, r, x) = a\a\_ v ..a\ 



,t 



(i + r + x) 2 (2r) 2 + £ ^ (^-j) 1 



—2— J 1 



It is not difficult to verify that the expression given in Eq. (fT2j) indeed solves Eq. (jllj) if one notes the recursion 
relation (fl"3"|) and the following properties 



K e (t, r, x) = 0, r 2 (d t + $.)X>(t, r, r - t) = - 



(2£ + l)! 1 
J\ (2r) e + 1 ' 
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of the integral kernel. Of course, the expression given in Eq. (fT2f is not the unique solution to Eq. (JTTJ ; one can add 
an arbitrary homogeneous solution to it. However, the solution is uniquely characterized by the following conditions: 

lim u^y, (const. + r, r) — (a^ vanishes at future null infinity), 
lim u^y, (const. — r,r) — ( u< e]ny vanishes at past null infinity), 

lim u^y,^, const.) = (^Hy vanishes at future time-like infinity). 

Notice that for £ > 1, the third requirement is necessary in order to exclude homogeneous solutions of the form 
*(t, r) — (c±t + c )r~ e (where cq and c\ are some constants), which vanish at both future and past null infinity. 
Summarizing, outgoing wave solutions have the form 

u lm Xt,r) = ty y^-^^ [^2(r - t) + 2Mc 3 U^\r - *)] + ^«£*"^>(t,r) +o(^£)\ 

(15) 

where the coefficients Cj are given in Eq. (| 1 3 j) and u (^ cfesca ' ter ) j g gi ven m Eq (|14p . A systematic derivation which 
includes the correction terms in M/R of arbitrarily high order is given in (l7j . 

Now let us analyze how much spurious reflection is introduced from these outgoing wave solutions when the condition 
Bl is imposed at the boundary surface. If for the moment we neglect correction terms arising from the backscatter as 
well as terms which are quadratic or higher order in M/R, then by the same arguments as in the previous subsection 
we conclude that Bl is perfectly absorbing for outgoing waves with angular momentum number I smaller or equal 
than L. Therefore, Bl automatically takes care of the curvature correction terms uf^ rv \ If effects from backscatter 
are taken into account, Bl is not perfectly absorbing for waves with £ < L anymore, however, in this case spurious 
reflections off the boundary surface are very small. In order to quantify this statement, consider for instance monopolar 
scalar radiation (i = 0, a = 0) for which 



oo 

/ \ TT , \ „„, f U 00 (x)dx ^/2M\ 

r-t 



(16) 



Since 



Uoojxjdx ( 2M\ 2M f U aQ (r-t + 2ry)dy ( 2M\' 
(d t + d r )uooXt,r)=-8M I {t + r + x)3 +0[- T ) =-—J + {-^) 

t o 



the boundary condition Bq is not perfectly absorbing unless M — 0. As a consequence, the solution of Eq. ([9j) consists 
of a superposition of an in- and an outgoing wave: 

woo(f,r) = u Q/>(t,r) + uoo\(t,r), 

where the ingoing wave has the form Uoa\(t,r) = Voo(r + t) + 0(2M/R) with Vbo a smooth function. Consider 
monochromatic waves of the form 

U m {r-t) = e k{t ~ r \ Vw{r + t)= 1 e-^ r+t \ (17) 

where k is a given wave number and 7 is an amplitude reflection coefficient. For the following, we assume that 
< k <gC M _1 . If kM is of the order of unity or larger, powers of kR which are multiplied by (2M/R) 2 might be 
comparable in size or larger than terms of the form M/R times unity, and in this case the (2M/i?) 2 -correction terms 
in Uooy(t,r) might actually be larger than the (2M/i?)-correction term. Imposing the boundary condition Bq at 
r = R, the ansatz (JT7J) yields 

Me 2lkR r e 2lkRy dy f2M\ 



7 R ikR J (1 + y) 3 + ° V R J 



i) 



It can be shown that the integral decays as (kR) 1 for large kR. Therefore, under the assumption that M <C k 1 <C R 
we find that I7I decays as (M/ R)(kR)~ 2 . Using similar arguments it can be shown that the boundary condition Bl 



6 



yields a reflection coefficient that decays as (M/R)(kR) ( L + 2 ) or faster for monochromatic waves satisfying I < L 
and M < k^ 1 < R. 

Finally, we remark that it is, in principle, possible to improve the boundary conditions Bl in order to take into 
account the first order correction terms in M/R of the backscatter. However, by the very nature of the backscatter, 
such boundary conditions cannot be local anymore. As an example, consider again monopolar scalar radiation where 
outgoing solutions have the form (fl6|) . Then, the boundary condition (l8l | 

t/2R 

fa , a\ u d\ , 2M f uoo(t~2Ry,R)dy 

(d t + d r )u Q0 (t, + J (1 + y) 3 = ®' ^ ' 

o 

with the boundary data 

oc 

nu . 2M f u 00 (0 ) R-t + 2Ry)dy 
G{t) -~lv J (TTW 3 ' 

t/2R 

is perfectly absorbing up to (and including) order 2M/R. Notice that the integral on the left-hand side of Eq. (fT8|) 
only involves the past portion { (r, R) : < r < t} of the boundary which is available from the past history of a 
Cauchy evolution starting at t = 0. The boundary data G involves an integral over the initial data at t = over 
the region r > R exterior to the computational domain and takes care of the backscatter that occurred in the past 
t < 0. If the initial data is compactly supported in the interval (0, R) this integral is zero and can be discarded. It 
was shown in [l8| that the boundary condition {T5J is stable in the sense that it admits an energy estimate. This 
construction can be repeated for waves with arbitrary angular momentum number I. 

A different method for constructing absorbing boundary conditions for linearized gravitational waves propagating 
on a Schwarzschild background has recently been presented in [2(| 0] . This method is based on fast converging 
series expansions of an exact nonlocal boundary condition and takes into account arbitrarily high correction terms in 
M/R of the Schwarzschild metric. However, there is no advantage to obtaining a boundary condition which takes 
into account the exact form of the Schwarzschild metric beyond the order of M/R in the construction of boundary 
conditions for wave propagation on a asymptotically flat curved background. The reason for this is that a generic, 
asymptotically flat background only agrees with the Schwarzschild metric up to order M/R. If second order effects are 
to be taken into account, quadratic terms in M/R and linear terms in J/R 2 (where J is the total angular momentum 
of the background) from the background metric must be considered. 



III. ABSORBING BOUNDARY CONDITIONS FOR EINSTEIN'S FIELD EQUATIONS 



The construction of absorbing outer boundary conditions in General Relativity is much more difficult than for the 
wave equation on a fixed background discussed in the previous section. At least three additional complications arise. 
First, in the Cauchy problem of General Relativity, Einstein's field equations split into a set of evolution equations 
and a set of constraints. If the spatial time slices are infinite or compact without boundaries it can be shown via 
the use of Bianchi's identities that any smooth enough solution of the evolution equations with constraint satisfying 
initial data automatically satisfies the constraints everywhere and at all times. However, if the time slices possess 
a nonempty boundary, this statement holds only if constraint-preserving boundary conditions are specified. The 
second complication is due to the fact that gravitational waves do not propagate on a fixed background but deform 
the spacetime metric as they evolve. As a consequence, it is not clear how to "fix" the boundary geometrically. It 
would be nice if one could specify the boundary conditions in such a way that the embedding of the boundary surface 
in the resulting spacetime is independent of the coordinate choice for which the evolution is performed. Otherwise, 
two evolutions using different coordinates might obtain different portions of spacetime even if they both start with 
the same initial slice and data. Finally, the third complication stems from the nonlinear nature of Einstein's field 
equations. In particular, the superposition principle for wave propagation does not hold so it is much harder to 
superpose an outgoing and and ingoing wave as was done in the previous section in order to quantify the amount of 
spurious reflection. 



A. Constraint-preserving boundary conditions 



The construction of constraint-preserving boundary conditions is probably the best understood and most studied 
issue of all the complications listed above: see Refs. @ HE 0, El, H 0, 0, 0, 
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IH, IM S3, SH H 0, for analytic studies and Refs. [H, [H, H HE [13, IH, H S3, H H S [HI, HI, [H, [H for 
numerical studies. The basic idea in constructing constraint-preserving boundary conditions is to derive the constraint 
propagation system, describing the propagation of constraint violations, and to impose boundary conditions for this 
system which ensure that zero is the only solution with trivial initial data. As an example, consider Einstein's field 
equations in harmonic coordinates, 

T a := U g x a = 0, 

where O g denotes the d'Alembertian operator with respect to the metric g ab . In these coordinates, Einstein's vacuum 
equations reduce to a set of ten coupled, quasilinear wave equations of the form 

g cd d c d d g a b = F ab {g)[dg, dg], 

subject to the constraint T a — 0, where F a i,(g)[dg, dg] depends quadratically on the first derivatives of the metric 
fields g a b- Using the twice contracted Bianchi identities, one finds that as a consequence of the evolution equations, 
r c satisfies a wave equation on its own, 

g cd d c d d T a = L a {g,dg,d 2 g)[T] 

where L a (g,dg,d 2 g)\T] is a linear first order differential operator in T a with coefficients depending on g ab and its 
partial derivatives up to second order. Therefore, the initial conditions r a | f=0 = 0, dtT a \ t=0 = insure that T a = 
on the domain of dependence of the initial slice. If the initial slice is a Cauchy slice, this implies that r a = everywhere 
on the spacetime, but if time-like boundaries are present, one needs to impose additional conditions at the boundary 
in order to guarantee T a = everywhere. There are different ways of assuring that T a = is the only solution with 
trivial initial data r a | 4=0 = 0, dtT a \ t=Q = 0. The simplest way is to specify Dirichlet data [4l[ 

However, other conditions, such as Sommerfeld-type conditions [3!j E3, HH, [Hj] or higher-order absorbing boundary 
conditions (56| are also possible. 

B. Gauge-controlling boundary conditions 

In the above example of Einstein's equations in harmonic coordinates, where the evolution equations have the form 
of ten coupled wave equations, one needs ten boundary conditions. As illustrated above, four conditions are needed 
in order to insure constraint propagation. Since there are two gravitational degrees of freedom, one expects that two 
boundary conditions are needed in order to control gravitational radiation. The remaining four boundary conditions 
are related to the residual freedom in choosing harmonic coordinates and fixing the geometry of the boundary surface. 
Preliminary ideas about how to specify such "gauge" controlling boundary conditions are given in [54] . 56] , but it is 
not clear yet how these conditions can be used in order to fix (at least part of) the geometry of the boundary surface. 
An exception are the boundary conditions constructed in [23[ for a tetrad formulation of Einstein's vacuum equations, 
which specify the mean curvature of the boundary surface as embedded in spacetime to be an arbitrary constant. 

Perhaps even more important than fixing the geometry of the outer boundary is the ability to specify a unique radial 
coordinate r and a unique outward radial null vector l a at each point of the boundary. Such quantities are needed 
for the generalization of the hierarchy of boundary conditions Bl defined in Eq. ©, where b — r 2 l a V a - A recent 
proposal for constructing r and l a based on the assumption that near the boundary, spacetime can be represented 
as Schwarzschild plus a small perturbation thereof is given in 18]. However, it is not yet completely clear how to 
identify the Schwarzschild background in this proposal. 

C. Absorbing wave boundary conditions 

Once constraint- and gauge-controlling boundary conditions have been specified, the next step is to construct 
boundary conditions which control the physical degrees of freedom by minimizing the amount of spurious gravitational 
reflection off the boundary surface. If the boundary is placed far from the strong field region, the field equations can 
be linearized about a weakly curved spacetime near the outer boundary. As discussed in Sect. Ill C[ it is sufficient to 
consider a Schwarzschild background provided R^> M where R is the radius of the outer boundary and M the total 
(Arnowitt-Deser-Misner) mass of the system. Therefore, near the outer boundary, it is safe to assume that spacetime 
can be written as Schwarzschild plus a small perturbation thereof. 
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Line ar p erturbations on a Schwarzschild background can be described by the Regge-Wheeler-Zerilli formalism 
[57L l58l. 153]. By performing a decomposition of the metric perturbations into spherical tensor harmonics, one obtains 
in this formalism two families of master equations for gauge-invariant potentials describing even (+) and odd 
(— ) parity metric fluctuations with angular momentum numbers £ > 2 and |m| < I. In particular, the metric 
perturbations in the Regge- Wheeler gauge can be reconstructed from the potentials < £ > ^ without solving additional 
differential equations [59l |. To first order in M/R, the master equations for have the form ([9]), where = 
—4 + 3/[(£ — l)(£ + 2)] in the even-parity case and a^' — —4 in the odd-parity case. Therefore, approximate 

= 0, where 



outgoing solutions have the form fT5]) and in principle, the boundary conditions Bl- fr i+1< ^m 



b = r 2 {dt + d r ), can be applied to the gauge- invariant quantities ■ However, the relation between and 
the metric perturbations is nonlocal in the sense that it depends on the angular momentum number £. Therefore, 
applying Bl in this way results in a nonlocal boundary conditions and its implementation requires a decomposition 
into spherical harmonics at the boundary. An alternative way is to first compute the (linearized) Weyl scalar 
from which is also a gauge-invariant quantity and to formulate the boundary condition on ^o- If (t ~ r^r) 

denote outgoing Eddington-Finkelstein coordinates for the Schwarzschild background, and = C a bcdl a m b l c m d is 
constructed from the Weyl tensor C a b c d and the null vectors l a d a — d t + d r , m a d a — (v^) -1 (<9,j + i sin -1 $ d v ) the 
boundary conditions [H, H3| 1 



-R 



s(±) 



Cr 



d t [r 2 (d t 



dr)} 



L-l 



r=R 



= 0, L = l,2,3,. 



(19) 



are perfectly absorbing for linearized gravitational waves with angular momentum number £ < L to first order in 
M/R if backscatter is neglected. Here, the time derivative operator dt in front of the operator inside the square 
brackets is introduced in order to allow for a static contribution to ^o- Notice that for L = 1 this condition just 
freezes V^o to its initial value. This freezin g-^ boundary condition has been given before in formulations of the 
IB VP of Einstein's equations [23L l36l |40L I42l.l47j . (Actually, the formulations in Refs. [23|, [36|, [4(| also consider more 
general boundary conditions which allow to couple ^4 to ^0 but it is unclear if this coupling is useful for reducing 
spurious reflections.) In this sense, the hierarchy Cl of boundary conditions improves the freezing- 'I' one. Reflection 
coefficients due to spurious reflections from backscatter for quadrupolar waves are computed in [181 ]. It is found that 
for a spherical outer boundary and quadrupolar gravitational radiation, Ci reduces spurious reflections by a factor 
of (15M/2R)(kR)~ 2 compared to the freezing-^o condition C\ when kR > 1. In 18], a new boundary condition 
T>2 similar to the the condition (fT5)) is derived which takes into account first order correction terms in M/R of the 
backscatter for quadrupolar linear waves. It would be interesting to generalize this analysis to take into account 
second order effects. However, in this case, it is not sufficient to consider perturbations of a Schwarzschild background 
since quadratic effects in M/R and linear effects in J/R 2 (where J is the total angular momentum of the spacetime) 
from the full metric should also be included. 



D. Well posedness results 



Once constraint-preserving absorbing boundary conditions with the desired properties have been specified, the next 
step is to prove the well posedness of the resulting IBVP. That is, one has to prove that for given initial data uo in 
an appropriate function space there exists a unique solution u(t) of the evolution equations in a time interval [0,T] 
which satisfies the constructed boundary conditions and such that u(0) = uo- Furthermore, one needs to show that 
u(t) depends continuously on the initial and boundary data in the sense that if — > uq is a sequence of initial data 
converging to uq then u^ n \t) converges to u(t) for each t G [0, T}. This property is important for the convergence of 
a numerical approximation since in this case the initial data always contains errors. Finally, one has to check that if 
Mo satisfies the constraint equation, so does u(t) for each t e [0, T]. 

A well posed IBVP for Einstein's vacuum equations was presented in Ref. [23j. This work, which is based on 
a tetrad formulation, recasts the evolution equations into first order symmetric hyperbolic quasilinear form with 
maximally dissipative boundary conditions [f5l], [62j for which (local in time) well posedness is guaranteed [63[ . There 
has been considerable effort to obtain well posed formulations for the more commonly used metric formulations of 
gravity. Partial results using similar mathematical techniques as in 23] were obtained in [24t |25| . |26I . l33l |34j, l38l . |43| . l44j . 



1 The normalization of the null vector l a differs from the one chosen in [Ti| by a factor of \f2/N , where N = 1 — 2M/r. This explains 
the absence of the factor N~ 1 in Eq. 1191 . 
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However, most of these works are either restricted to the linearized equations or to reflecting (and not absorbing) 
boundary conditions. For results on coupled hyperbolic-elliptic linear problems with constraint-preserving boundary 
conditions based on semigroup techniques see [35], H3] • 

A different technique for showing the well posedness of the IBVP is based on the frozen coefficient principle 
where one freezes the coefficients of the evolution and boundary operators. In this way, the problem is simplified 
to a linear, constant coefficient problem on the half-space which can be solved explicitly by using a Fourier-Laplace 
transformation [64] . This method yields a simple algebraic condition (the determinant condition) which is necessary 
for the well posedness of the IBVP. Work based on verifying the determinant condition for the Einstein case is given 
in [22], [U l33L l34l [H, HH ■ Sufficient conditions for the well posedness of the frozen coefficient problem were developed 
by Kreiss [65|. Kreiss' theorem provides a stronger form of the determinant condition whose satisfaction leads to 
well posedness if the evolution system is strictly hyperbolic. One of the key results in [65[ is the construction of a 
smooth symmetrizer for the problem for which well posedness can be shown via an energy estimate in the frequency 
domain. Using the theory of pseudo-differential operators it is expected that the verification of Kreiss' condition also 
leads to well posedness for quasilinear problems, like Einstein's field equations. Work based on the verification of 
Kreiss' condition in the Einstein case is given in [H, \^ but since in that case the evolution system is not strictly 
hyperbolic it is not clear if these results imply well posedness. Recently, Kreiss and Winicour [41[ introduced a new 
pseudo-differential first order reduction of the wave equation which leads to a strictly hyperbolic system. Using this 
reduction they were able to prove well posedness of the IBVP for Einstein's field equations in harmonic coordinates 
in the frozen coefficient approximation. This work is generalized to higher-order absorbing boundary conditions in 
[56| . For an alternative proof of the results in (4lj which does not require a pseudo-differential first order reduction 




E. Applications to numerical relativity 



For applications of constraint-preserving boundary conditions to numerical relativity, see Refs. [451 [46[ for simula- 
tions of self-gravitating scalar fields in spherical symmetry, Ref. [47( for the simulation of ID collidi ng grav itational 
plane waves, Ref. [5l[ for evolutions of Brill waves in axisymmetry, Refs. 0, [H, 128^ 36^ 37l . l39l [HoL [H, l53l |54| for tests 
in three spatial dimensions, Ref. [55| for binary black hole simulations and Refs. 67T6q for the simulation of bubble 
spacetimes in five-dimensional theories of gravity. In particular, Refs. (36l . I37I ] implement a first order symmetric hy- 
perbolic formulation of Einstein's vacuum equations with the freezing-\l/o boundary condition, and Refs. [1^, HH, [55| 
also freeze \l/o at the boundary but use the harmonic formulation of the field equations. In (54[, constraint-preserving 
freezing-^ boundary conditions are tested for the case of a perturbed Schwarzschild black hole and compared to 
other types of boundary conditions proposed in the literature. It is found that the version of constraint-preserving 
freezing-^o boundary conditions in [54} performs better than all alternate boundary treatments tested. It should be 
interesting to numerically implement the boundary conditions Cl, L > 2, given in Eq. (|19|) which are refinements of 
the freezing-^o boundary condition. 



IV. CONCLUSIONS 



Formulating absorbing outer boundary conditions for the numerical solution of Einstein's field equations involves five 
steps: i) The construction of constraint-preserving boundary conditions which make sure that no constraint- violating 
modes enter the computational domain, ii) finding boundary conditions that geometrically control the evolution of the 
boundary surface, iii) finding conditions that minimize the amount of spurious reflection of gravitational radiation off 
the boundary, iv) proving well posedness of the resulting initial-boundary value problem (IBVP) and v) discretizing 
the problem. 

As discussed in this article, there has been a lot of effort in carrying out step i) which is, by now, well- understood. 
In contrast to this, step ii) needs further work. Regarding step iii), a promising approach for minimizing spurious 
reflections is the hierarchy Cl of local boundary conditions on the Weyl scalar ^0 presented in [l8| . They have the 
property of being perfectly absorbing including curvature corrections (but neglecting backscatter) to order M/R for 
all multipoles of gravitational radiation up to L, where M is the Arnowitt-Deser-Misner mass of the spacetime and 
R a typical radius of the boundary surface. However, their precise formulation requires a radial coordinate and an 
outward radial null vector field at the boundary whose unambiguous definition is an open problem and could benefit 
from progress in step ii). Regarding step iv), a complete proof of the well posedness of the IBVP has been given in [23| 
for a tetrad formulation of Einstein's field equations. Recently, proofs for well posedness have also been given for the 
frozen coefficient limit of Einstein's equations in harmonic coordinates (4ll . |5(| and it is expected that these results can 
be generalized to the full Einstein equations based on the theory of pseudo-differential operators. Finally, in step v), 
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promising results have been achieved in the numerical implementation of the harmonic formulation with constraint- 
preserving absorbing boundary conditions [H, [H, [55| . It is expected that the boundary condition C 2 proposed in [r| , 
which is perfectly absorbing for quadrupolar linearized gravitational radiation, will improve these results. Finally, it 
would be interesting to develop similar boundary conditions for the Baumgarte-Shapiro-Shibata-Nakamura [o9l |70T ] 
formulation of Einstein's field equations which is often used in numerical relativity. For partial results along these 
lines, see @H EL 
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